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^ ■ Abstract 

, The ellipsoidally symmetric Buda-Lund hydrodynamic model describes 

CO ' naturally the transverse momentum and the pseudorapidity dependence 

' of the elliptic flow in Au+Au collisions at y^ijvjv = 130 and 200 GeV. The 

result confirms the indication of quark deconfinement in Au+Au collisions 
at RHIC, obtained from Buda-Lund hydro model fits to combined spectra 
and HBT radii of BRAHMS, PHOBOS, PHENIX and STAR. 



Introduction. PHENIX, PHOBOS and STAR experiments at RHIC pro- 
duced a wealth of information on the asymmetry of the particle spectra with 
Q ' respect to the reaction plane P EJ 13 ^ ISJ |B] , characterized by the second har- 

^ , monic moment of the transverse momentum distribution, denoted by W2- It is 

C ' measured as a function of the transverse mass and particle type at mid-rapidity 

as well as a function of the pseudo-rapidity rj = 0.5 log(|^p^). 
►v>( \ The PHOBOS collaboration found 3 , that V2{ri) is a strongly decreasing 

; I ' function of |?7|, which implies that the concept of boost-invariance, suggested by 

\ Bjorkcn in ref. 7 , cannot be applied to characterize the hadronic final state of 

Au+Au collisions at RHIC. 

We summarize here a successful attempt to describe the pseudo-rapidity 
dependence of the elliptic flow V2{'ri) at RHIC, for more details see ref. Our 
tool is the Buda-Lund hydrodynamic model OEII, which we extended in ref. |H] 
from axial to ellipsoidal symmetry. 

Buda-Lund hydro for ellipsoidal expansions. Based on the success of the 
Buda-Lund hydro model to describe Au+Au colhsions at RHIC jElESI, Pb+Pb 
colhsions at CERN SPS and /i + p reactions at CERN SPS [HIISl, we 
describe the emission function in the core-halo picture, and assume that the 
core evolves in a hydrodynamical manner: 

o / ^ j4 9 P^d'^Y.^{x) 

Sc{x,p)d X f- — , (1) 

{2TTy B(x,p) + Sq 
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where g is the degeneracy factor [g = I for identified pseudoscalar mesons, g = 2 
for identified spin=l/2 baryons), and p^d^Y,^{x) is a generahzed Cooper-Frye 
term, describing the flux of particles through a distribution of layers of freeze- 
out hypersurfaces, B{x,p) is the (inverse) Boltzmann phase-space distribution, 
and the term Sq is determined by quantum statistics, Sq = 0, —1, and -1-1 for 
Boltzmann, Bose-Einstein and Fermi-Dirac distributions, respectively. 

For a hydrodynamically expanding system, the (inverse) Boltzmann phase- 
space distribution is 

Bi-^P)=^-P[-Yi^-T^)- (2) 

We will utilize some ansatz for the shape of the flow four- velocity, Ui,{x), chem- 
ical potential, /i(a;), and temperature, T(x) distributions. Their form is deter- 
mined with the help of recently found exact solutions of hydrodynamics, both 
in the relativistic [T7] and in the non-relativistic cases El 120] ■ 

The generalized Cooper-Frye prefactor is determined from the assumption 
that the freeze-out happens, with probability iJ(r)dT, at a hypersurface char- 
acterized by r = const and that the proper-time measures the time elapsed 
in a fluid element that moves together with the fluid, dr — M''(x)(ia;^. We 
parameterize this hypersurface with the coordinates {rx,ry,rz) and find that 
(i'^I]^(x|T) = u'^{x)d'^x/u'^{x). Using dtT\r = u'^{x) we find that in this case the 
generalized Cooper-Frye prefactor is 

p^d^i:^{x) =p^u^{x)H{T)d^x, (3) 

This finding generalizes a result of ref . "21" from the case of a spherically symmet- 
ric Hubble flow to anisotropic, direction dependent Hubble flow distributions. 

From the analysis of CERN SPS and RHIC data [HI [HI 021 , we find that 
the proper-time distribution in heavy ion collisions is rather narrow, and H{t) 
can be well approximated with a Gaussian representation of the Dirac-delta 
distribution, 

with At <C tq. 

We specify a fully scale invariant, relativistic form, which reproduces known 
non-relativistic hydrodynamic solutions too, in the limit when the expansion 
is non-relativistic. Both in the relativistic and the non-relativistic cases, the 
ellipsoidally symmetric, self- similarly expanding hydrodynamical solutions can 
be formulated in a simple manner, using a scaling variable s and a corresponding 
four- velocity distribution u^, that satisfy 

u'^d^.s = 0, (5) 

which means that s is a good scaling variable if its co-moving derivative van- 
ishes [13 nH]. 
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It is convenient to introduce the dimensionless, generalized space-time ra- 
pidity variables {r]x, rjy, rjz), defined by the identification of 

sinhrjx^r^^, (6) 

similar equations hold for y and z. The characteristic sizes (for example, the 
lengths of the major axis of the expanding ellipsoid) are {X, Y, Z) that depend 
on proper-time t and their derivatives with respect to proper-time are denoted 
by {X, Y, Z). Eq. © is satisfied by the choice of 

coshr7r — 1 coshn,, — 1 cosh 77, — 1 

= (7, sinh77:r, sinhriy, sinh rjz), (8) 

and from here on {Xf,Yf,Zf) = {X{to),Y{to), Z{Toj) = {Xi,X2,X3), assum- 
ing that the rate of expansion is constant in the narrow proper-time interval of 
the freeze-out process. The above form has the desired non-relativistic limit, 



^2 

V 



2X] 2Yf 2Zy 



(9) 



where again {Xf,Yf,Zf) = (X(ro), y(ro), Z(ro)) = [X^^X^.X^). From now 
on, we drop subscript /. The normalization condition of u'^{x)Ufj_{x) = 1 yields 
the value of 7. For the fugacity distribution we assume a shape, that leads to 
Gaussian profile in the non-relativistic limit, 

Pf\-^-s, (10) 
T(x) To 

corresponding to the solution discussed in refs. |18[ 1191 1^ . We assume that 
the temperature may depend on the position as well as on proper-time. We 
characterize the inverse temperature distribution similarly to the shape used in 
the axially symmetric model of refs. [51 110|. and discussed in the exact hydro 
solutions of refs. ^1 E| j 



T{x) To V J \ Te 2A 

where Tq, T, and Te are the temperatures of the center, and the surface at 
the mean freeze-out time tq, while Te corresponds to the temperature of the 
center after most of the particle emission is over (cooling due to evaporation 
and expansion). Sudden emission implies Te = Tq and At — > 0. 

The observables can be calculated analytically from the Buda-Lund hydro 
model, using a saddle-point approximation in the integration. This approx- 
imation is exact both in the Gaussian and the non-relativistic limit, and if 
p'^u^/T ^ 1 at the point of maximal emittivity. 
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The results are summarized in Figs. 1 and 2. We find that a small asymmetry 
between the two transverse Hubble constants gives a natural explanation of the 
transverse momentum dependence of V2 ■ The parameters are taken from Buda- 
Lund hydro model fits in refs. |12[ \'2'2\ . where the axially symmetric version of 
the model was utilized, although the mean transverse flow is reduced due to 
the less central nature of collisions studied in this study, which is a summary 
of ref. [Sj. Note that we did not yet fine-tune the ellipsoidal Buda-Lund hydro 
model to describe these V2 data, instead we searched the parameter set by hand. 
So at the moment we are not yet ready to report the best fit parameters and the 
error bars on the extracted parameter values. However, as indicated by Figs. 
1 and 2, even this fit by eye method is successful in reproducing the data on 
elliptic flow at RHIC. 

First we tuned the model to describe the pt dependent elliptic flow of identi- 
fied particles at midrapidity, as shown in Fig. 1. Then we calculated the value 
of the transverse momentum integrated V2 {r] = 0) and found [H) , that this value 
is below the published PHOBOS data point at mid-rapidity. We attribute the 
difference of 0.02 to a non-flow contribution The PHOBOS collaboration 

pointed out the possible existence of such a non-flow contribution in their data 
in ref. as they did not utilize the fourth order cumulant measure of V2- 

We note that our presently best choice of parameter set correspond to a 
high, Tq > Tc — 170 MeV central temperature, with a cold surface temperature 
of Ts « 105 MeV, see Figs. 1 and 2. 

Summary and conclusions. We have generalized the Buda-Lund hydro model 
to the case of ellipsoidally symmetric expanding flreballs. We kept the param- 
eters as determined from fits to the single particle spectra and the two-particle 
Bose- Einstein correlation functions (HBT radii) |12l 122] , and interpreted them 
as angular averages over the direction of the reaction plane. Then we observed 
that a small splitting between the expansion rates parallel and transverse to 
the direction of the impact parameter, as well as a small tilt of the particle 
emitting source is sufficient to describe simultaneously the transverse momen- 
tum dependence of the collective flow of identifled particles as well as the 
pseudorapidity dependence of the collective flow [2113 at RHIC. 

The results confirm the indication for quark deconfinement at RHIC found 
in refs. |12[I22| . based on the observation, that some of the particles are emitted 
from a region with higher than the critical temperature, T > Tc = 170 MeV. The 
size of this volume is about 750 fm'^, corresponding to 1/8-th of the total volume 
measured on the t = tq main freeze-out hypersurface fS] . At the same time, this 
analysis indicates that the surface temperature is rather cold, Tg « 105 MeV, so 
approximately 7/8 of the particles are emitted from a rather cold hadron gas. 
So the picture is similar to a fireball, which is heated from inside. 
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Figure 1: Buda-Lund fit to the V2{pt) data 
Here we see the fit to the PHENIX V2 (pt) data of identified particles [T|. The parameter 
set is: To = 210 MeV, X = 0.57, Y = 0.45, Z = 2.4, Ts = 105 MeV, ro = 7 fm/c, 
■0 = 0.09, Xf = 8.6 fm, Yf = 10.5 fm, Zf = 17.5 fm, fio,-^ = 70 MeV, fio,K = 210 MeV 
and /io,p = 315 MeV, and the masses are taken as their physical value. 




Figure 2: Buda-Lund fit to the V2{ri) data 
This image shows the fit to the 130 GeV Au+Au and 200 GeV Au+Au V2{rj) data 
of PHOBOS EHH, with the ellipsoidal generalization of the Buda-Lund hydro model. 
Here we used the same parameter set as at fig. Q with pion mass and chemical 
potential, and a constant non-fiow parameter of 0.02. 
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